Instability of isolated triplet excitations on the Shastry-Sutherland lattice (SSL) 
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Configurations of singlets and triplets on the SSL have been proposed in the literature as varia- 
tional ground states of the Shastry-Sutherland model at fixed magnetization M. We prove, that iso- 
lated triplet excitations on the SSL are unstable if the coupling a falls below a critical value ac ~ 2.0. 
The instability should be visible in the compound SrCu2 (603)2 where a coupling a* — 1.48 is re- 
alized. 



PACS numbers: 75.10.-b, 75.10.Jm 
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INTRODUCTION 

The discovery of plateaus in the magnetization curve 
M = M{B) of the compound SrCu2(B03)2 § |,i at ra- 
tional values of the magnetization M /Ms = 1/3, 1/4, 1/8 
{Ms = 1/2) has led to intensive investigations of the 
ground state properties of the Shastry-Sutherland model 
[01 with Hamiltonian 



H = S(x)S(y)-Ka ^ S(x)S(y) 

(x,y> ((x,y)) 



(1) 



The spin- 1/2 couplings extend over nearest ((x, y)) and 
specific next- nearest (((x, y))) neighbor couplings, which 
define the SSL, shown in Fig. ^ 



where plateaus are expected or found - several ordering 
patterns have been proposed and discussed in the litera- 
ture. ||l|, H, H, 0, |, ^ They all have in common that 
in the sector with total spin S*, Nt = S triplets and 
Ns — N/2 — Nt singlets are distributed over the SSL 
with N sites. Typical examples of these singlet-triplet 
configurations arc shown in Rcf ^. In contrast to the 
pure singlet configuration (Fig. |l|) the singlet-triplet con- 
figurations are not eigenstates of the Hamiltonian p). 
Moreover, it has been shown in our previous paper llG] 
that the formation of triplets on the SSL is disadvan- 
tegeous if the coupling a falls below a certain critical 
value ac{M), depending on M. In particular, the singlet- 
triplet configurations proposed in Fig. 8 of Ref. |^ cannot 
be considered as adequate variational states for M = 1/8, 
since ac{M = 1/8) = 2.3 > a*. As an alternative, we 
proposed in Ref. ^ (Fig. 4) a "monomer-dimer configu- 
ration" built up from monomers (i.e. isolated spin-up 
states) and dimers (i.e. singlets) on nearest and next- 
nearest neighbor bonds. 

It is the purpose of this Letter, to prove that the in- 
stability of isolated triplet excitations on the SSL is a 
quite general feature for a < 2 and low magnetizations 
M < 1/8. 



SINGLET-TRIPLET VERSUS 
MONOMER-DIMER CONFIGURATIONS 



FIG. 1: The couplings in the Shastry-Sutherland model. 
Nearest and next-nearest neighbor couplings are represented 
by dotted and solid lines, respectively. The diagonal bonds 
define the Shastry-Sutherland lattice (SSL). 

In the absence of an external field (A/ = 0) the ground 
state can be represented by a product of singlet states on 
the SSL if a exceeds a critical value etc, which has been 
found to be ac — 1.43 in Ref. ||, ^. This value is very 
close to the coupling a* ~ 1.48, which seems to be real- 
ized in the compound SrCu2 (603)2. In order to describe 
the magnetic properties of the compound - in particu- 
lar for fixed magnetizations M = 1/4, 1/6, 1/8, 1/12, 1/16 



The variational states \K, v) which are constructed 
from the monomer-dimer configurations generalize the 
familiar singlet-triplet configurations on the SSL in the 
following sense: In the sector with total spin S we have 
11 = 23 monomers - i.e. isolated spin-up states |xj-|-) 
j = 1, . . . , J/ at sites Xi . . . x^. They form the ferromag- 
netic cluster {v). All remaining sites x, y, . . . are occupied 
with {N — v)/2 singlets (dimers) [x, y]. They form the 
antiferromagnetic cluster {K). Each monomer and dimer 
configuration defines a variational state: 

V 

\K,v) = nix,+) n t^'y]- (2) 



The expectation value of Hamiltonian (|^) between these 
states turns out to be [[lOl 



{K,iy\H\K, v) 



where n[°\k) and N\''{y) are the numbers of nearest 
neighbor singlets on {K) and monomer pairs on {y)^ re- 



(1) 



+ aN!f\K) 



{u)+aN!^\v) 



(3) 



r(i) 



spectively. N^\K) and N2'''{v) denote the correspond- 
ing numbers on the next-nearest neighbors according to 
Fig. [|. The singlet-triplet configurations on the SSL are 
special monomer-dimer configurations with: 



(0) 



0, N, 



(0) 



(1) 



0, N. 



(1) 



(4) 



i.e. we have a maximum number of dimcrs (singlets) 
and monomer pairs (triplets) on the SSL. Due to the 
geometry of the SSL it is not possible to maximize the 
number of sing lets N^f^ and to minimize the number of 
triplets n!^^ at the same time. Here begins the problem 
with triplet excitations on the SSL: For smaller values 
of a [a < ac{M)\, the break up of the triplets into well 
separated monomer pairs is favored energetically. For 
M < 1/8 - where the triplets on the SSL are well isolated 
(cf. Figs. 8,9,10 in Ref. ^) - we can study the break up 
of each triplet and the change of the local environment 
in the monomer-dimer configuration directly, as is shown 
in Fig. ||(a),(b). The configuration {Ka^Va) in Fig. ||(a) 
shows a triplet at sites x and y/ embedded in a sea of 
singlets on the SSL. They are represented by the solid 
lines in Fig. H(a). The triplet is broken up in a pair 
of monomers at sites x and y separated by a "knight's 
move" [cf. Fig. |2|(b)]. Comparing the energies of the 
configurations (_ft'a,;/a) and {Kf^i^b) we have to decrease 



Af2°'' and N2''' and to increase N^"' by one unit: 
AE = {Kb,J^b\H\Kb,iyb) ~ {Ka,iya\H\Ka,jya) 



(0) 



2^-2 



(5) 



The energy difference ^ is negative for a < 3/2. There- 
fore, the break up of the triplet is favored in this regime. 
We will see in the next section, that this phenomenon 
occurs already at larger values of a, if wc improve our 
variational ansatz. 



MAGNETIC ORDER AND THE FROZEN 
MONOMER APPROXIMATION 

The monomer-dimer configurations are not eigenstates 
of the Hamiltonian (|l|). However, it might happen that 
the ground state is governed by a specific distribution 
(i/) of "frozen" monomers. They define the magnetic 



(a) 



(b) 




(c) 
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FIG. 2: (a) An isolated triplet on the SSL. The singlets (solid 
lines) and the dotted couplings define the antiferromagnetic 
cluster (Kaji'a)- (b) Break up of the triplet into isolated 
monomers at sites x and y and the antiferromagnetic cluster 
{Kb, Vb)- (c) The antiferromagnetic 8-site cluster, which leads 
to the energy difference (bl). 



ordering of the ground state. The variational ansatz (g) 
will be improved if we substitute the product of dimers 



n [x,y] ^ ^{K) 



(6) 



by the ground state "^{K) of the antiferromagnetic clus- 
ter Hamiltonian H{K): 



H{K)'^{K) = E{K)^{K). 



(7) 



H{K) contains all nearest and next-nearest neighbor 
couplings on the antiferromagnetic cluster K. In 
Fig. ^(a),(b) they are represented by the solid singlet 
lines and the dotted lines connecting them. 

The variational ansatz (|^) yields an upper bound on 
the ground state energy Eq{M = v/2N) < E{K,v) in 
the sector with magnetization M, where 

E{K,j.) = \ (N^^\j.)+aN!,'\v))+E{K). (8) 

Concerning the triplet-singlet configuration in Fig. ^a), 
the ground state "^{Ka) of the antiferromagnetic cluster 



3 



Hamiltonian H{Ka) is again given by the product of sin- 
glets on the SSL. Note that each singlet (e.g. [1, 2]) is ac- 
companied by two nearest neighbor couphngs S(l)S(z), 
S(2)S(2:) such that the total spin operator S(l) + S(2) 
acts on the singlet [1, 2]. For this reason the ground state 
energy cannot be lowered in the case (i^a, i^a)- 

The situation is different in configurations with iso- 
lated monomers. The ground state ^'(i^f,) of the an- 
tiferromagnetic cluster is not identical with the dimcr 
product shown in Fig. [2(b)]. If we apply here the cou- 
plings S(2)S(3), S(4)S(5) and S(4)S(7) onto the nearest 
neighbor singlet [3,4] a new state is created. Therefore, 
the ground state energy E{Ki,) of the antiferromagnetic 
cluster Kb will be lower than the expectation value of the 
corresponding monomer-dimer configuration in Fig. ^(b). 
This will shift the critical value ac where the triplet ex- 
citation becomes unstable to larger values ac > 3/2. In- 
deed, we can derive an improved lower bound for ac, if 
we substitute the dimers [1, 2] [3, 4] [5, 6] [7, 8] by an anti- 
ferromagnetic cluster with 8 sites 1,...,8 as shown in 
Fig. §(c). The ground state energy £'§(«) of the 8-site 
cluster is lower than the energy of the corresponding 
dimer configuration in Fig. ^(b) by an amount 

Asia) = -^{3a+l)-Es{a); (9) 

As{a) is shown in Fig. ^. The straight line AE — {a ~ 
3/2)/2 represents the energy difference (||). Therefore, 
the intersection point 

ac = 1.95 (10) 

yields the critical value, where the instability of the 
triplet sets in. Enlarging the antiferromagnetic surround- 
ing of the frozen monomer pair in Fig. y(c) leads to an 
increase of the critical value (|lo|). However, the effect is 
small: a calculation with a 16-site cluster - containing 
the 8-site cluster - yields ac = 1.97. 

The changes in the magnetic order from an isolated 
triplet on the SSL to a pair of monomers separated by 
a "knight's move" can be tested by measuring the spin- 
spin structure factors, which are obtained from spin-spin 
correlators 5'3(x)53(y). 

The expectation values of the spin-spin correlators can 
be read off from the monomer-dimer configurations: 

{K,y\Ss{^)Ss{y)\K,,) = i{5W(x,y)-5(")(x,y)} 

(11) 

where 6'^^'>{x,y) = 1 and ^^''^(x, y) — 1 only if (x,y) co- 
incides with nearest or next-nearest neighbor bonds on 
{K, v) occupied with a monomer pair or a dimer, respec- 
tively. In all other cases (5^^^ (x, y) — 0. Note in particu- 
lar, that all correlators with one site x on the ferromag- 
netic cluster V and the other site y on the antiferromag- 
netic cluster K vanish. 



In the frozen monomer approximation with an anti- 
ferromagnetic cluster wavefunction 5* (if), (|l|) can be 
extended to 

{^(K),v\S:,{x)S^{y)\^{K),v) = (12) 
i<5(i)(x,y) + (vI/(/0|53(x)53(y)|vI/(if))5(x,y e K) . 

Note that ( pT| ) and (|l^) differ in the correlations 
53(x)53(y), x, y G K on the antiferromagnetic cluster 
K. 
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FIG. 3: The energy differences @M and The inter- 

section points define the critical coupling ac where the break 
up of the triplet sets in. 



COEXISTENCE OF STABLE AND UNSTABLE 
TRIPLET EXCITATIONS ON THE SSL FOR 

M = 1/6 

The instability discussed so far only concerns "iso- 
lated" triplets on the SSL. The density of triplets 
Nt/N = M is controlled by the magnetization AI and 
raises the question, what happens in the cases M = 1/4 
and M = 1/6. A typical singlet-triplet configuration for 
M = 1/4 on the SSL is given in Fig. 6(a) of Ref. |. 
Breaking up a triplet in these configurations produces 
additional monomer pairs, which increases the energy, 
i.e. a high density of triplets prevents their instabil- 
ity. In the case M = 1/6, the triplets on the SSL form 
stripes as is shown in Fig. ^(a). Interesting enough, 
each second stripe of the triplets can break up as is 
demonstrated in Fig. ^(b). Comparing the energies of 
the configurations {Ka,Va) and {Ki,,Vb), we find again, 
that the energy difference (||) per triplet changes sign for 
adM = 1/6) = 3/2. 

Note, that the isolated monomers along one stripe in 
Fig. ^(b) are accompanied with a quasi-onedimensional 
antiferromagnetic cluster K. We have computed the 
ground state energy of the 24 site clusters E24{a,M = 



4 



1/4), which surrounds the 4 pairs of monomers arising 
from the break up of 4 triplets on the SSL, shown in 
Fig. |^(b). The energy difference (per triplet) 



A24(a,l/6) 



?a + 4) -S24(a,l/6) W13) 



is plotted in Fig. ||. It meets the straight line (a — 3/2) /2 
at 



adM = 1/6) = 1.985. 



(14) 



This improved lower bound for ac{M — 1/6) differs from 
the lower bound we found in Ref. There we did not 
realize, that a mixture of triplets on the SSL and isolated 
monomers - as shown in Fig. ^(b) - lowers the energy 
more efficiently. 
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FIG. 4: (a) Singlet-triplet configuration (Ka,Va) on the SSL 
for M = 1/6. (b) Break up of triplets along a stripe with an 
antiferromagnetic 24-site cluster {Kb,iyb)- 



DISCUSSION AND CONCLUSIONS 

In this Letter, we have investigated the stability of 
triplet excitations on the Shastry-Sutherland lattice. Our 



results can be summarized as follows: 

(i) Isolated triplets turned out to be unstable if the 
ratio of next-nearest to nearest neighbor couplings is 
smaller than a < ac — 2.0. This implies that the mag- 
netic properties of SrCu2 (603)2 - where the couphng 
ratio is a* ~ 1.48 - cannot be understood in terms of 
isolated triplet excitations on the SSL, if the magneti- 
zation M is small enough (M < 1/8). This statement 
holds in particular for the triplet ordering patterns for 
M = 1/8, 1/12, 1/16 in Figs. 8,9,10 of Ref. |. 

(ii) For Af = 1/6 [Fig. 7(a) in ReL | the triplets form 
stripes on the SSL. It turns out that the triplets on each 
second stripe break up into isolated monomers if a < 
ac(M = 1/6) = 1.985. 

(iii) The triplets on the SSL in the configuration for 
M — 1/4 [Fig. 6 (a) in Ref. ^ appears to be stable for 
a = a* . This demonstrates, that triplet excitations on 
the SSL are stabilized if the triplet density is high enough. 

Breaking up the isolated triplet into two monomers 
(spin-up states) at sites x and y separated by a "knight's 
move" [Fig. §(c)] lowers the energy of the new configura- 
tion in two ways: 

1. There is a "classical" effect, which results from the 
rearrangement of singlets and triplets on (Kajiya) 
[in Fig. ^(a)] into the monomcr-dimcr configuration 
{K.iyh) [in Fig. 1(b)] . 

2. There is a "quantum" effect produced by the 
antiferromagnetic cluster, which surrounds the 
monomers at sites x and y [Fig. |^(c)] . 

We are indebted to M. Karbach for a critical reading 
of the manuscript. 
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